Hydrodynamic effects in interacting Fermi electron jets 
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Abstract 

We theoretically study hydrodynamic phenomena originating from electron-electron collisions 
in a two-dimensional Fermi system. We demonstrate that an electron beam sweeping past an 
aperture creates a pumping effect, attracting carriers from this aperture. This pumping effect 
originates from the specific electric potential distribution induced by the injected electrons. In the 
regions nearby the main stream of injected electrons, a positive potential is induced by the injected 
electrons. Thus, the normally repulsive Coulomb interaction leads to an attractive force in the 
Fermi system. This quantum pumping mechanism in a Fermi system differs qualitatively from the 
Bernoulli pumping effect in classical liquids. We also discuss possible experimental realizations. 
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Electron-electron (e-e) scattering in a degenerate electron gas forms one of the central 
concepts of Fermi-liquid theory^* 2 . Electron states with energy slightly above the Fermi 
surface are not stationary, since they decay accompanied by the creation of electron-hole 
pairs in the Fermi seai. The total momentum in the e-e scattering process is conserved, and 
hence the e-e collisions do not in the main affect the conductivity. At the same time, the e-e 
collisions play an important role in experiments which involve the phase coherence time of 
the electron^, and in physical effects implicating hydro dynamic o 5 ' 6 ' 7 > 8 . At low temperatures, 
the e-e collision processes slow down due to Fermi statistics ^ 10 ' 11 , leading, in the ideal two 
dimensional electron system (2DES), to the inverse e-e lifetime rapidly decreasing upon 
lowering the temperature T^°< 

L = <ML lbi .EjL +bi _±_ + (1) 

r ee hE F 1 k B T a*k F J ' v ; 

where h, E F , and v F are the Planck's constant, the Fermi energy, and the Fermi veloc- 
ity, respectively; denotes the effective Bohr radius. To render hydrodynamic effects in 
2DESs observable, e-e collisions should occur more often than impurity scattering events. 
Mathematically, this condition is obvious: l ee < l p , where l ee = v F r ee is the electron mean 
free path related to e-e scattering and l p is the momentum relaxation mean free path. The 
above condition is satisfied in typical high- mobility GaAs/GaAlAs heterostructures, in the 
temperature range T ~ 5 — 35 K (Fig. la)^. Therefore, the hydrodynamic effects should be 
accessible in experimental low-temperature studies. Here we describe a new hydrodynamic 
phenomenon in a low-temperature electron plasma which occurs when an electron beam is 
injected from a narrow aperture into a 2DES (Fig. lb). We show that the injected beam 
creates an unusual potential pattern which can result in a hydrodynamic pumping effect. 
This pumping effect exists in mesoscopic structures where the injected electrons experience 
just few e-e collisions. In contrast to the Bernoulli effect in classical liquids, the predicted 
quantum pumping behaves linearly in the excitation, and can only be observed in a Fermi 
system at low temperatures. 

Another potential application of the predicted hydrodynamic effect concerns measure- 
ments of the electron phase coherence time, r^. The low-temperature behavior of in 
nanostructures is a long standing problem. According to the standard theory, e-e scattering 
becomes very weak at low temperatures and therefore should tend to infinity at zero 
temperature. However, the experiments demonstrate that saturates to a finite value at 
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low temperatures^. Magnetic-impurity scattering has been implicated amongst others^. 
Yet, e-e scattering itself is a nontrivial problem in mesoscopic structures, where an electron 
experiences scattering by boundaries as well as impurities. For example, it is accepted that 
in diffusive metals, the formula l/r ee oc T 2 should be strongly modified due to disorder—. 
Measurements of the hydrodynamic effect predicted in this work can directly reveal the 
contribution of e-e scattering to the phase coherence time in mesoscopic ballistic structures 
at low temperatures. 

To describe hydrodynamic effects in a stationary electron beam (Fig. lb) we employ the 
linearized Boltzmann equation: 

v^ + eE^ = J (/l ), (2) 
or op 

where /i(r, p) is the distribution function of the non-equilibrium electrons, r = (x,y) and 
p are the two-dimensional coordinate and momentum, respectively; /o( e ) is the equilibrium 
Fermi function, v = p/m is the electron velocity, and e and m are the energy and the 
effective mass of the electron, respectively; E is the in-plane electric field, e = — \e\ is the 
electron charge, and J is the collision integral. At low temperature, the non-equilibrium 
electrons have energies close to the Fermi surface and the function f\ can be written as 
/i(r;e, #) = — ^-F(r;9), where 9 is the angle between the velocity and the direction +x, 
with < 9 < 2n. Our choice of (x, y) coordinates is indicated in Fig. lb. 

We treat e-e collisions in the relaxation-time approximation. Due to the complexity of 
the collision integral incorporating the long-range Coulomb interaction, theoretical studies 
based on the long-range Coulomb collision integral usually include only the first e-e collision 
event^iia. In contrast, the relaxation-time approach allows us to obtain analytical results for 
an infinite number of collisions and to describe fundamental hydrodynamic effects appear- 
ing in the multi-collision regime. Within the relaxation-time approximation, the collision 
integral describing scattering events in the vicinity of the Fermi surface, takes the form: 



F - F - 2 cos #cos 9F - 2 sin #sin 9F 
J(F) = (3) 



where F = / F(6)d9/2<K and r ee is the e-e scattering time. Note that the e-e collision 
integral conserves both the number of particles and the momentum. In addition, the 
collision integral (jSJ) assumes an ideal system without disorder and is valid if / ee < l p . 
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To solve the Boltzmann equation in the half-plane x > we need to impose the boundary 
condition (BC) at the line x = 0. Outside the aperture, the BC describes the elastic collision 
of an electron with an ideal border, F(0, y; 9) — F(0, y; tt — 9) =0. In order to express the 
injection process, we now introduce the injection function g(y, 9), nonzero only in the angle 
intervals < 9 < n/2 and (3/2)7r < 9 < 2ir. It is convenient to consider the BC separately 
in two intervals, —it/2 < 9 < tt/2 and 7r/2 < 9 < (3/2)n. In the first angular interval, the 
outgoing flux of electrons vfF(9) should equal the sum of the incoming and the injected 
fluxes, vpF^ii — 9) + vpg{9). A similar argument is applied to the second interval. In both 
intervals, the BC condenses to a single formula: 

F(0, y- 9) - F(0, y- n - 9) = g(y; 9) - g(y; n-9). (4) 

Importantly, the BC (@J) leaves some freedom in choosing the function F(0,y; 9). A general 
solution of eq. EJcan be written as F(0, y \ 9) = Wi(y; 9) + w 2 (y; 9), where the function Wi(9) 
is determined by the injection g(y;9), whereas the function W2{9) is arbitrary and satisfies 
the elastic-collision condition u>2(y; 9) — W2(y; n — 9) = 0. 

The BCs in this problem are given by the equation 0] at the line x = and by the 
condition F(x — > oo) — > at infinity. To solve this problem in a rather convenient way, 
we now consider a symmetric problem, involving the entire two-dimensional plane and two 
symmetric current sources (two jets) in the vicinity of r = 0. Including the two current 
sources, we write the linearized Boltzmann equation as follows: 

dF 

v— - eEv = J(F) + G(y; 6)6 (x) + G(y; n - 9)6 (x), (5) 
or 

where the function G(y; 9) describes the injection. By exploiting the symmetry of the prob- 
lem and integrating eq. El over x in the vicinity of x — 0, we determine that the solution of 
eq. El satisfies the necessary BC (HJ) with G{9) = vp cos 9g{9). Equation (JSJ) can be solved 
using Fourier transformation, 



• i i7 F k ~ F k - 2 cos #cos 9F k - 2 sin 6>sin 9F k 

ivkF k - eE fc v = + G ky (9) + G ky (ir-9), (6) 

where Fk(9) = J^°° <i 2 re _jkr F(r; 9). The electric field E in eq. |H1 originates from the non- 
equilibrium electrons. We now assume that the heterostructure containing the 2DES is 
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covered by a top metallic gate, with the distance between the 2DES and the top gate 
d smaller than any lateral size of the system. This condition does not unduly restrict 
experimental situations. Imposing this condition ensures that the potential in the plane of 
2D plasma is proportional to the non-equilibrium 2DES density 5n(r). Thus, we can write 
0(r) = (47red/e s )Sn(r), where e s denotes the dielectric constant of the semiconductor—^. 
After Fourier transformation, we obtain = — ik0fc = — faulted/ e s )8rik. Simultaneously, 
the electron density is expressed through the function F, 5n k = F^Z^d, where D 2 d = m/7rh 2 
represents the 2D density of states. By manipulating eq. |H1 and integrating over angles, we 



obtain a closed system of equations for cos 6^, and sin^-FV The 2DES density can be 
determined by the functions 

= = /lffi(W0-l)+/2 m 

k (l + 4d/o5)(l-A;oWo)' [) 

wm = -* w= r G f^J M (8 ) 

a/ k 2 + k\ Jo nv F J ik cos{6 - a) + k 2nv F 

where tana = k y /k x and k = {vft^)" 1 = l/l ee . We also derive analytical expressions for 
the current densities, j x (k) = VFD 2 DCOs9F k and j y (k) = VFD 2 Dsin9F k . 
For the injection function, we choose 

G(y,») = G m „e-^5(^_M, (9) 

with 0(0) = 1 if 9 > and otherwise. The parameters L and 9$ describe the spatial width 
of the aperture and the width of the angular distribution of injected electrons, respectively. 
The parameter G max is proportional to the injected current density in the middle of the 
aperture. Using the injection function, G max can be related to the total injected current, 
by I — \e\D2DL/(2y/7r)G max . Typically, the resistance of the injecting aperture exceeds 
the resistance of the leads and thus the potential drop across the aperture, AV, dominates 
the potential applied to the system. Using the Landauer-Biittiker formalism, we can write 
I = (2e 2 /h)NAV—, where iV represents the number of conducting modes in the aperture. 
Although quantization of the conductance is not necessary for the observation of the hy- 
drodynamic effect (N can be large), the Landauer-Biittiker formalism forms a convenient 
framework to link the experimental variables, L, N, AV and I. By combining the above 
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equations, we obtain: G max = 4y/7r\e\N/(hLD2D)AV. In addition, it was assumed that the 
injected electrons possess energies in the vicinity of the Fermi surface. 

Numerical results for the non-equilibrium density distribution follow from the reverse 
transformation, Sn = D2D J exp(ik.r)Fkdh/(27r) 2 . Figure 2 shows that Sn is large and posi- 
tive inside the ballistic beam, as expected. But Sn is depressed to negative values in regions 
nearby the main stream of injected electrons. The negative Sn corresponds to depletion. 
In Fig. 3, we schematically show the streamlines associated with injection. Interestingly, in 
the regions adjacent to the main stream, the currents flow towards the beam. Qualitatively, 
we can understand this behavior in terms of e-e scattering: the injected electrons create an 
effective pressure and scatter the background Fermi-sea electrons toward the right. Numeri- 
cal calculations show that this happens not in the main stream, where the density of excess 
electrons is high, but rather on the sides of the main stream. The currents flowing toward 
the main stream tend to compensate for the lack of electrons. 

The electron depletion in the vicinity of the injected beam can lead to carrier pumping 
toward this beam. The potential induced by the non-equilibrium electrons can be expressed 
as 0(r) = (47red/e s )Sn(r), reflecting the fact that a local net charge will lead to a local 
potential of the same sign. If a detector aperture is situated in the region of electron 
depletion (Sn < 0), the detector will experience a positive potential (Fig. 3). In voltage 
measurement mode, the detector can be considered a closed reservoir with net zero current. 
A positive potential denotes a lack of electrons in the detector reservoir, the result of a net 
pumping of electrons toward the injected beam. We can also consider the pumping in terms 
of currents. Suppose that we turn on the injected current I at t — 0. In the region of the 
detector window, the pumping effect extracts electrons from the detector. After some time, 
the system reaches steady state, and the net current through the detector window vanishes. 
Hence a counter current of electrons into the detector must be generated, by a positive 
potential on the detector lead. If on the contrary current is allowed to flow, electrons will 
be pumped through the aperture in steady state. To achieve the steady state pumping, we 
can connect the detector contact to the Fermi sea in the right-hand side, including a resistor 
to maintain the current at a level sufficiently low for the latter to be regarded as a weak 
perturbation. In all above cases, the positive potential induced by the hydrodynamic effects 
in the vicinity of the injected beam, effectively result in a pumping phenomenon. 

According to eq. HJ the length l ee = VFT ee strongly depends on temperature, leading 
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to a temperature dependence of the pumping effect. In the ideal disorderless system, the 
pumping potential depends on the ratios L/l ee , x/l ee , and y/l ee (Fig- lb). It is revealing to 
analyze the temperature dependence of the induced electric potential at a fixed position with 
coordinates (xa,Ua) where a detector window can be located (Fig. 3). Calculations (Fig. 3) 
indicate that the potential induced by the injected beam at the point A with coordinates 
(xa = 2 fim, Da = —2 fim) is negative at low and high temperatures, while for intermediate 
temperatures, 3 — 25 K, the hydrodynamic effects render the potential positive. In this 
regime, l ee ~ xa ~ Va and electron-electron interaction leads to just a few scattering events 
in the vicinity of the detector. Therefore, the pumping effect occurs only in the regime of few 
scattering events. Our calculations show that the pumping voltage increases by decreasing 
the angle 8 . Hence, better experimental conditions to observe the pumping signal can 
be found in relatively wide apertures with flared openings, where the collimation effect is 
typically stronger 17 . In addition, since oc 5n oc G max oc AV oc J, the pumping voltage is 
a linear function of the applied voltage or injected current. 

In the case of one-dimensional injection (L — > oo), the function 5n(x) can be obtained 
from eq. El by setting Gk y oc 5{k y ). The calculated electron density clearly demonstrates the 
hydrodynamic effect: the Fermi sea depletes in the region x ~ l ee due to the effective pressure 
created by the injected beam (Fig. 4). This effect can also be applied toward a pumping 
mechanism involving a quasi-one- dimensional channel featuring an injecting barrier and a 
detector window, as illustrated in Fig. 4. 

It is interesting to compare our results with the hydrodynamics of classical liquids. Of 
course, the local non-linear equations of fluid mechanics differ qualitatively from the non- 
local kinetic Boltzmann equation. For example, the Bernoulli pumping effect, resulting from 
the spatial variation of local speed in the moving liquid, is quadratic in terms of fluid speed, 
hence non-linear. As result the Bernoulli pumping effect remains unchanged if we reverse 
the current. In the 2DES case discussed in this work, the pumping force is linear in terms of 
applied voltage or current and changes sign as we reverse the current flow. The linearity is 
a consequence of the electron Fermi statistics. Indeed, at low temperatures, we can consider 
the flow equally in terms of electrons or of holes near the Fermi surface. Mathematically, 
the described linear response originates naturally from the linearized Boltzmann equation. 
Physically, when the negative voltage is applied to the injecting contact, the injected beam 
consists of electrons which induce an effective pumping force for the electrons from a detector 
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(Fig. 3). When a positive voltage is applied, the injected beam can be treated as a beam of 
holes near the Fermi surface; these holes now induce an effective pumping force for the holes 
from the detector and therefore the pumping voltage reverses its sign, along with the applied 
voltage. The concurrent voltage reversals lead to a linear dependence at low excitations. 
The linear dependence stands in contrast to classical liquids where currents consist of like 
particles and the above Fermi-gas description cannot be applied. 

To conclude, we have shown that the Coulomb interaction in an injected ballistic beam 
results in attractive forces which can be exploited toward a pumping effect. The quantum 
pumping described by us is a peculiar property of a degenerate 2DES and is qualitatively 
different from the hydrodynamic effects in classical liquids. 

We thank Axel Lorke for helpful discussions. This work was supported by NSF (DMR- 
0094055) and by the Condensed Matter and Surface Science Program at Ohio University. 
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FIG. 1: 

(a) Momentum relaxation and electron-electron mean free paths, l p and I function of 
temperature. The parameter L is obtained from mobility measurements on a GaAs/AlGaAs 
heterostructure at the 2DES density 4.8* 10 11 cm~ 2 [13]. I ee is calculated using eq. 14 from ref. 10. 

(b) Geometry of the system with an injected electron beam, (c) Schematic of two symmetric jets 
in a two-dimensional plasma. 
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FIG. 2: Calculated non-equilibrium electron density as a function of the in-plane coordinates. 
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Schematic of a hydrodynamic pump in a mesoscopically patterned 2DES. The stream lines and 
density distribution are shown schematically. The detector serves as a probe for the potential near 
the main beam of electrons. Insert: The induced voltage at the point A with coordinates (xa, Ua) 
as a function of temperature; number of conducting modes in the QPC, N = 10; 2DES density = 
4.8 * 10 n cm~ 2 , d = 400 A, a* = 100 A, and e = 12.5. 
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FIG. 4: 

Calculated non-equilibrium electron density as a function of the x-coordinate in the one-dimensional 
system. Insert: Schematic of a hydrodynamic pump utilizing a one-dimensional channel with an 
injecting barrier and a detector aperture. 



13 



